MULTI-COMPONENT GAS ANALYZING METHOD USING FTIR 

Field of the Invention 

The present invention relates to a multi-component gas analyzing method using FTIR 
(Fourier Transform Infrared Spectroscopy), which executes a quantitative analysis on a plurality 
of components in a sample based upon an absorption spectrum obtained by FTIR. 

Description of Prior Art 

FIG. 1 is a schematic view showing the structure of a gas analyzer (hereinafter, referred 
to as FTIR gas analyzer) that uses a Fourier transform infrared spectroscopy conforming to 
FTIR. In FIG. 1, reference number 1 is an analyzing section, reference number 2 represents a 
data processing section for processing interferrograms, which is an output of the analyzing 
section 1 . 

The analyzing section 1 is constituted by an interference mechanism 7 consisting of an 
infrared light source 3 arranged to emit parallel infrared light rays, a beam splitter 4, a fixed 
mirror 5 and a movable mirror 6 that is allowed to make parallel movements, for example, in X- 
Y directions by a driving mechanism installed outside FIG. 1, a cell 8 in which a test sample, a 
reference sample, etc. are placed and to which the infrared light rays from the infrared light 
source 3 are directed through the interference mechanism 7, and a semiconductor detector 9, etc. 

Here, the data processing section 2, which is formed by, for example, a computer, is 
arranged to average the interferrogram, transform the averaged interferrogram into absorption 
spectrum by fast Fourier transformation and carry concentration calculations on target 
components based upon the absorption spectrum. 

In the FTIR gas analyzer arranged as described above, a plurality of components are 
quantitatively analyzed in the following manner. Reference samples or test samples are 
respectively flowed into the cell 8, and infrared rays from the infrared light source 3 are directed 
to the cell 8 so that interferrograms of the reference samples or the test samples are measured. 
These interferrograms are respectively Fourier-transformed in the data processing section 2 to 
obtain a power spectrum. Then, the ratio of the power spectrum of the test sample to the power 
spectrum of the reference sample is calculated, and an absorption spectrum is obtained by 
converting the ratio into the absorbance scale. Thereafter, a plurality of components contained in 
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the test sample are quantitatively analyzed based upon absorbances of a plurality of wave 
number points in the absorption spectrum. 

With respect to the method for quantitatively analyzing a plurality of components, for 
example, the applicant of the present application have proposed several methods in the 
5 specifications of the following Japanese Patents (JP Nos.: 2649667, 2741376, 2926277, and 
2926278). For example, JP No. 2649667 discloses a multi-component analyzing method in 
which the sum of relative absorbances, each representing a difference between a local peak value 
and a local valley value in each point of a plurality of wave number points in the absorption 
spectrum, is found. Based upon the sum, concentrations of a plurality of components are 
10 respectively obtained independently. 

In addition to the above-mentioned methods, with respect to the quantitative analysis on 
,^ multi-components, quantitative algorithms for multi-variable analysis, such as PLS (Partial Least 

O Square), PCR (Principal Component Regression) etc., may be used. In any of these methods, 

Ul 

Sj concentration operating data (calibration matrix) are preliminarily prepared based upon a 



03 



y A 5 calibration spectrum having a known concentration, and by applying the calibration matrix for a 
sample spectrum having an unknown concentration, concentrations of a plurality of components 
in the sample spectrum are calculated individually. 

However, in the case of measurements on engine exhaust gases by the use of the above- 
mentioned FTIR, it is generally known that if the generally used quantitative algorithm, such as 
10 PLS, PCR, etc., is applied, the readings of CO, C0 2 and NO become higher than the actual 
concentrations due to an influence of high concentration of H 2 0 commonly existing in a sample 
gas at the time of an actual measurement of the sample gas. This is not due to the interference 
that occur by an insufficient separation of overlapped spectra even by the use of the quantitative 
algorithm, but is due to the fact that there is a difference in the spectrum intensity itself, for 
25 example, as shown in FIG. 6, even in the case of the same component and the same 
concentration, depending on the base gas composition. More specifically, this is because CO in 
an actual sample gas which contains both N 2 and H 2 0, as it were (N 2 + H 2 0) base gas, shows a 
greater absorption spectrum than CO in the span gas containing only N 2 as a base component. In 
FIG. 6, reference numbers 10, 11 and 12 show the respective absorption spectra of CO, (CO + 
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H2O) and H2O when N2 gas is used as a base gas. Additionally, the same phenomenon can be 
observed even in the coexistence of O2, H2 which have no infrared-ray absorption. 

Although an actual mechanism of the above-mentioned phenomenon is not clearly 
understood, it is considered that quenching due to interactions between gas molecules might 
5 relate to this phenomenon as one of the reasons. FIG. 7 is a model illustrating changes in the 
amount of infrared absorption due to quenching, and this model indicates that changes in the 
amount of infrared absorption occur depending on the degree of probable collision between a gas 
component X and a base gas and the interaction between them. In other words, FIG. 7(A) shows 
a case in which both the possibility of collision and the interaction of a base gas A to the gas 
10 component X are small. In this case, hardly any effects are given to the equilibrium of the basic 
condition and excited condition of gas component X so that the existence of the base gas A 
hardly gives any effects on the amount of absorption of the gas component X. Moreover, FIG. 

P 

5 7(B) shows a case in which both the possibility of collision and the interaction of the base gas B 
l\ to the gas component X are great. In this case, since the equilibrium or the gas component X is 
Cfl5 shifted to the basic condition side, additional light absorption is more likely to occur. In other 
jf{ words, the existence of the base gas B increases the absorption intensity of the gas component X, 
® thereby providing a stronger absorption as compared with the case of the base gas component A. 
p In the example shown in FIG. 6, N2 represents components A, H 2 0 represents component B and 
CO represents component X. 

Moreover, another possible mechanism is "collision expansion" in which the line width 
Q of the absorption spectrum is widened due coexistence with another component. 

Summary Of Invention 

The present invention has been devised to solve the above-mentioned problems, and its 
objective is to provide a multi-component gas analyzing method using FTIR, which can correct 
25 for an influence of a coexistent gas due to an intensity change in the spectrum itself that occurs 
depending on the base gas composition so that it is possible to carry out a measurement with 
high precision. 

In order to achieve the above-mentioned objective in a multi-component gas analyzing 
method using FTIR of the present invention, a plurality of components in a sample are 
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quantitatively analyzed based upon an absorption spectrum obtained by FTIR. After having 
calculated multi-component concentrations from a mixed gas spectrum by using a quantitative 
algorithm, calculations are further carried out so as to correct for a change in spectrum due to a 
coexistent gas component with some target components. 

5 Upon carrying out the above-mentioned correcting calculations, if the coexistent gas 

component is measured by FTIR, the measured results can be used. In this case, there is no time 
delay among the measured values. Moreover, if the coexistent gas component is not measured 
by FTIR, the corresponding data may be read from an external analyzer and used, and in this 
case, time matching may be made in the CPU of the FTIR. 

10 In accordance with the multi-component gas analyzing method using FTIR of the present 

invention, it is possible to correct an error caused by an intensity change in the spectrum itself 
which can not be considered by a normal quantitative algorithm, and consequently, to improve 
C the precision of a gas analysis using FTIR. In particular, the present invention is effectively used 
Cj for correcting influences due to a difference in the base gas composition between an exhaust gas 
9 5 and a calibration gas, that is, in particular, a difference in H 2 0 concentration and a change in the 
y3 H 2 0 concentration in a sample gas being measured. 

JU Brief Description Of The Drawings 

FIG. 1 is a schematic view showing one sample of a device for carrying out the method 
of the present invention; 

Q>0 FIG. 2 are exemplary graphs showing the influence of coexistent of H 2 0 with respect to 

the CO 10% range, wherein FIG. 2(A) shows the case of fixed CO concentration and FIG. 2(B) 
shows the case of fixed H 2 0 concentration; 

FIG. 3(A) is a graph that shows a model of the influence on CO and H 2 0, and FIG. 3(B) 
is a graph that shows one example of a correction expression; 
25 FIG. 4 are graphs showing the relationship between the CO concentration (calculated 

values) and the error, wherein FIG. 4(A) shows a state before the correction and FIG. 4(B) shows 
a state after the correction; 

FIG. 5 are graphs showing variations in the influenced value in the CO 10% range, 
wherein FIG. 5(A) shows actual data and FIG. 5(B) shows data after having been standardized; 
30 FIG. 6 is a graph that explains conventional problems; and 
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FIG. 7 is a model illustrating changes in the amount of infrared absorption due to 
quenching. 

Detailed Description Of The Preferred Embodiments 

Referring to FIGS. 1-7, the following description will discuss embodiments of the present 
5 invention. As described above, the measurement of an actual sample gas is greatly influenced by 
H2O contained therein. Therefore, an explanation will be given of an exemplary of examination 
on correction expression correcting the influence of coexistent H2O which is a component that 
can be measured by FTIR on CO or other components after concentration calculations carried 
out by a multivariable analyzing method. Here, it is taken as a premise that influences caused by 
10 components other than H 2 0 are actually negligible. 

1. Examination on Correction Expression 

(1) Confirmation of Trend of Influenced Value 

FIG. 2 are exemplary graphs showing the influence of coexistent H2O with respect to the 
7fi CO 10% range (calculated by raw concentration values before linearizing), wherein FIG. 2(A) 
tfl5 shows the case for a fixed CO concentration and FIG. 2(B) shows the case for a fixed H 2 0 
S concentration. As shown by FIG. 2 for the fixed H2O concentration, the influence can be 
approximated by a linear equation (Y = aX). In some cases, depending on the components, the 
i influence may be better approximated by a quadratic equation (Y = aX + bX). 

if (2) Model of Influence and Correction Expression 

3 

20 FIG. 3 shows examples of a model of influence and a correction expression in CO and 

H2O. In this model, it is supposed that the peak intensity of H2O (nearly equal to calculated 
concentration value before linearizing) x (=W) and a change in the peak intensity of CO (a 
change in concentration calculated value before linear rise) y have a relationship represented by a 
linear equation (y = Sx). The inclination S of this linear equation changes in proportion to the 

25 peak intensity of CO as shown in FIG. 3(A). The relationship between the inclination S (=Y) 
and the peak intensity X (=C) of CO without influence of H2O is approximated by a quadratic 
equation (Y = aX 2 + bX), as shown in FIG. 3(B). Here, the CO concentration value R, obtained 
by the conventional PCR, is represented by: 

R = y + X- {(aX 2 + bX)Xx} + X 



in 
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By solving this equation with respect to X (=C), a relational expression for finding the 
peak intensity of CO in which the influence of H2O is corrected is obtained. 

In the above-mentioned method, in order to obtain correction coefficients a and b, data 
shown in FIG. 3(B), that is, actual measured values of the influence of coexistence at the time 
when only the CO concentration is changed with the H 2 0 concentration being virtually fixed 
(data at the time when the H2O concentration is changed are not necessary), are used. 

Then, the CO concentration (before linearizing) C after correction is represented as 
follows: 

whena = 0(Y = bX), 

C = R/(l+bW) 
when a * O (Y= aX 2 + bX), 
C = [ V {(1 + bW) 2 + 4aWR} - (1 + bW)]/2aW 
where R: CO concentration (before linearizing) prior to correction, 

W: H 2 0 concentration (before linearizing). 

(3) Example of Result of Correction 

FIG. 4 are graphs showing an example (CO 10% range) of correction made by the 
expression shown in FIG. 2, wherein FIG. 4(A) shows data before the correction and Fig. 4(B) 
shows data after the correction. As a result of this correction method, the influenced value, 
which is approximately 12% at maximum (full scale, FS), is maintained within a range of ±1 % 
or slightly more. 

2. Confirmation of Difference in Devices 

It is necessary to decide how to take actual data for determining coefficients for each of 
devices based upon variation of the absolute value and trend of the influence of coexistent H 2 0. 
Therefore, three devices were used so as to confirm the variations and the possibility of 
simplification for the correction method. 

(1) Confirmation of Trend of Influenced Values 
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FIG. 5 are graphs showing variations in influenced values in the CO 10 % range, wherein 
FIG. 5(A) shows the relationship between the CO peak intensity and errors in actual data (CO 
10%) and FIG. 5(B) shows the relationship between the CO peak intensity and errors in a 
standardized state (CO 10%). As shown in FIG. 5, although there are considerable gaps in terms 
of absolute values, there is a common trend. It is considered that one of the reasons for 
variations is a difference in wave number points used in calculation for PCR in each of the 
devices. 

Moreover, the following Table 1 shows variations in influence in which other 
components were included. Table 1 indicates that CO, CO2, NO and N2O are subjected to great 
influences, and with respect to these, since the width of variations is great, the correction 
coefficient needs to be determined for each device. 
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[Table 1] 



Component 


Range 


Influence 


(error, 


FS%) 


Average 


Maximum 












(FS%) 


Deviation 














n7<S%\ 






A 


B 


C 






CO 


200ppm 


6.3 




f.J 


0.4 






1 uuuppm 


O.J 


7 3 






-1 1 






0 0 


8 1 
0 . 1 


1 1 9 

1 1 ..7 


9 6 


2.2 




9% 


Q Q 


8 O 


13 7 


10.5 


3.2 






1 9 7 


10 3 


18 0 


n 7 


4.3 




1 % 
1 /o 


1 0 ^ 


1 1 0 


8 0 


.7.0 


-1 8 






9 0 


7.3 


7.4 


7.9 


1.1 




20% 


11.7 


11.9 


3.6 


9.0 


-5.5 


IN \J 




3 7 


3.1 


5.0 


3.9 


1.1 




1 OOOnnm 


6 7 


6.0 


9.5 


7.4 


2.1 




5000Dnm 


8.7 


9.0 


11.4 


9.7 


1.7 


N 2 0 


200Dt»m 


5.1 


5.2 


7.2 


5.8 


1.4 


ism. 


500nnm 


5.1 


3.4 


2.3 


3.6 


1.5 


so 2 


200ppm 


-.24 


0.9 


-4.4 


-1.9 


2.9 


CH30H 


500ppm 


3.6 


3.2 


1.9 


2.9 


-1.0 




2000rrom 


2.4 


1.5 


0.7 


1.5 


-0.9 


CH 4 


500ppm 


1.9 


1.2 


1.9 


1.7 


-0.4 


C 2 H 4 


500ppm 


2.8 


1.1 


-3.1 


0.3 


-3.3 


C3H6 


200ppm 


-0.9 


1.0 


-2.3 


-0.7 


1.7 


1,3-C4H6 


200ppm 


2.1 


1.2 


-1.6 


0.6 


-2.2 


iso-C 4 H 8 


200ppm 


-0.4 


1.5 


-3.5 


-0.8 


-2.7 



(2) Simplification of Correction 

With respect to the measurement results of the three devices A to C shown in Table 1, an 
explanation will be given of the results of corrections made on components, CO, CO2 and NO. 
In these cases, the following three correction methods are suggested: 

(a) A quadratic equation (a * O in FIG. 3) is used as a correction coefficient, and this is 
determined for each of the devices. 

(b) A quadratic equation (a * O in FIG. 3) is used, and based on the correction coefficient 
of one of the three devices A, B, C, the inclination is adjusted for each of the devices. 
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(c) A linear equation (a = O in FIG. 3) is used as a correction coefficient, and this is 
determined for each of the devices. 

The following description will discuss the respective correction methods and the results 
thereof: 

5 (a) In the case where the correction coefficient (quadratic equation) is determined for 

each device. 

Five data points (1/5 division point to 5/5 division point) of Wet (N 2 + H 2 0 base) span 
spectra were sampled for each of the components of CO, C02, and NO. Based upon the data 
points, the correction coefficient was determined for each device. 

10 The results thereof are shown in Table 2. 

[Table 2] 



m 


Component 


Range 


Influenced 


Value 


(FS%) 








A 


B 


C 


CO 


CO 


200ppm 


0.1 


-1.2 


0.3 


i. ™ 




lOOOppm 


0.2 


0.4 


0.3 


ft* 




5000ppm 


0.2 


0.1 


0.1 






2% 


0.1 


-0.1 


0.2 






10% 


0.1 


0.1 


0.3 




C0 2 


1% 


0.1 


0.2 


0.3 


ru 




5% 


0.1 


0.2 


0.2 


yi 




20% 


0.1 


0.1 


-0.2 


Q 


NO 


200ppm 


0.3 


-0.4 


0.4 


O 




lOOOppm 


-0.1 


0.2 


0.2 






5000ppm 


0.2 


0.3 


0.2 



(b) In the case where only the inclination is adjusted by using the same coefficient 
(quadratic equation). 

15 Among data of the respective devices used in Table 2, the correction coefficient 

(quadratic equation) for device A is revised and used based upon a single data point (4/5 division 
point). The results from the corrections of the remaining four data points are evaluated. The 
results thereof are shown in Table 3. Here, the results related to device A are the same as those 
in Table 2. 
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[Table 3] 



Component 


Range 


Influenced 


Value 


(FS%) 






A 


B 


c 


CO 


200nom 

a* \y \*r Ly 1/111 


0.1 


-2.4 


-0.8 




lOOOppm 


0.2 


-1.2 


-0.5 




5000ppm 


0.2 


-0.2 


-0.3 




2% 


0.1 


0.3 


-0.6 




10% 


0.1 


0.2 


0.9 


C0 2 


1% 


0.1 


-0.2 


-1.0 




5% 


0.1 


-0.3 


-0.2 




20% 


0.1 


0.3 


-0.5 


NO 


200ppm 


0.3 


-0.6 


-1.0 




lOOOppm 


-0.1 


-0.3 


-0.5 




5000ppm 


0.2 


1.0 


-0.4 



(c) In the case where the linear equation is used to correct each of the devices. 

The correction coefficient (linear equation) was determined by using only a single data 
point (4/5 division point) among data of the respective devices used in Table 2. The results from 
correcting the remaining four points are evaluated. The results thereof are shown in Table 4. 



[Table 4] 



Component 


Range 


Influenced 


Value 


(FS%) 






A 


B 


C 


CO 


200ppm 


-0.6 


-1.9 


-0.8 




lOOOppm 


-0.5 


-1.6 


-1.0 




5000ppm 


-0.4 


-0.4 


-0.6 




2% 


-0.6 


-0.4 


-1.2 




10% 


-0.7 


0.5 


-1.4 


C0 2 


1% 


1.3 


-1.5 


-1.1 




5% 


0.2 


-0.4 


-0.3 




20% 


0.3 


-0.3 


-0.6 


NO 


200ppm 


-0.6 


-0.8 


-0.9 




lOOOppm 


-0.6 


-0.8 


-0.7 




5000ppm 


-0.6 


-0.6 


-0.8 
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Tables 2 to 4 show that, even in the case of the linear equation (c), the results of the 
correction are maintained within a width of approximately ±2%. Consequently, it is possible to 
make corrections even by sampling only a single data point of actual data for each 
component/range. 

As described above in the present invention, in a multi-coefficient gas analyzing method 
using FTIR in which a plurality of components in a sample are quantitatively analyzed based 
upon an absorption spectrum obtained by FTIR, after having calculated multi-component 
concentrations, from a mixed gas spectrum by using a multi-variable analyzing method, etc., 
calculations are further carried out so as to correct a change in spectrum due to a coexistent gas 
component with respect to some target components. Therefore, it is possible to correct an error 
caused by an intensity change in the spectrum itself which can not be corrected by a normal 
quantitative algorithm and, consequently, to improve the precision of a gas analysis using FTIR. 

In particular, for example, in the case of H 2 0 where the correction is made by using 
measured values obtained in FTIR, the method is advantageous because the time delay among 
measured values for the respective components are negligible. Moreover, with respect to 
components such as O2 and H 2 that are components which FTIR fails to measure in principle, 
signals taken from an external analyzer can be time-matched in the CPU of the FTIR and then 
used in the correction. Thus, it is possible to obtain the same effects. 

The multi-component gas analyzing method of the present invention is effectively used 
for correcting the influence caused by a difference in water concentration between exhaust gases 
and a calibration gas and a change in H2O concentration in a sample gas being measured. The 
method also makes it possible to correct influences caused by O2 and H2. 
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